Nonclassical correlation in a multipartite quantum system: two measures and evaluation 
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There is a commonly recognized paradigm in which a multipartite quantum system described by a density 
matrix having no product eigenbasis is considered to possess nonclassical correlation. Supporting this paradigm, 
we define two entropic measures of nonclassical correlation of a multipartite quantum system. One is defined as 
the minimum uncertainty about a joint system after we collect outcomes of particular local measurements. The 
other is defined by taking the maximum over all local systems about the minimum distance between a genuine 
set and a mimic set of eigenvalues of a reduced density matrix of a local system. The latter measure is based on 
an artificial game to create mimic eigenvalues of a reduced density matrix of a local system from eigenvalues of a 
density matrix of a global system. Numerical computation of these measures for several examples is performed. 
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I. INTRODUCTION 

There has been a longstanding discussion on the definition 
of quantumness in a quantum state of a multipartite system. 
One definition is, of course, entanglement which is consid- 
ered equivalent to inseparability according to the separability 
paradigm |Q], IU, . The separability paradigm suggests, as is 
well-known in this field, a classification of density matrices of 
a system consisting of subsystems 1, . . . , m. Separable den- 
sity matrices are those of the form 

k 

with positive weights Wk (%2k Wk = 1) anc ^ density matrices 
p]} of subsystems. Inseparable density matrices are those that 
cannot be represented in this form. A system (consisting of 
remote subsystems) represented by a separable density matrix 
is regarded as a classically correlated system because local op- 
erations and classical communications (LOCC; see, e.g., [4]) 
can create it from scratch: It can be prepared remotely when 
distant persons (Alice, Bob,. . ., Marry) receive instructions 
from a common source (Clare). 

A bipartite system is a typical system to investigate. Bipar- 
tite separable density matrices are those of the form 

k 

with positive weights Wk (%2k w k — !)■ Bipartite inseparable 
density matrices are those that cannot be represented in this 
form. One supporting evidence for the paradigm is that, bi- 
partite system represented by a separable density matrix does 
not violate Bell's inequality 12, HI]. 
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Detection methods of inseparability have opened a large re- 
search field, many of which are based on the Peres-Horodecki 
test Hfl using positive but not completely positive linear 
maps. 

It is still a challenging issue to find classes of density matri- 
ces possessing nonlocal nature, other than the class of insep- 
arable density matrices. Bennett et al. (6D discussed a certain 
nonlocality about locally nonmeasurable separable states. Ol- 
livier and Zurek [7] later introduced a measure called quan- 
tum discord defined as a discrepancy of two expressions of 
a mutual information that should be equivalent to each other 
in a classical information theory. Another branch of study on 
quantumness was started by Oppenheim and the Horodecki 
family [8]; this was extensively studied by a group consisting 
of the Horodecki family and other authors 190 . They intro- 
duced a protocol called closed LOCC (CLOCC). This pro- 
tocol allows only local unitary operations, attaching ancillas 
in separable pure states, and operations to send subsystems 
through a complete dephasing channel. They also defined a 
measure of quantumness named quantum deficit as a discrep- 
ancy between the information that can be localized by apply- 
ing CLOCC operations and the total information of the sys- 
tem. The present work is in the stream of these studies on 
quantumness in correlation. 

One way to evoke a discussion on the validity of the sep- 
arability paradigm is to look at the persistent question why a 
pseudo-pure state 

ft» = Pl^X^I + (i - P)l/d (3) 

with an entangled pure state and d the dimension of the 
Hilbert space, is often regarded as a classically correlated state 
for small probability p because of its separability proved by 
Braunstein et al. [10]. The state p ps can be regarded as a state 
possessing quantumness in correlation if we choose another 
paradigm than the separability paradigm. It is thus a rather 
conceptual question and is in relation to the following discus- 
sion. 

Suppose that we have a system consisting of two sub- 
systems (local systems) and cannot eliminate a local and/or 
global superposition by local unitary operations. A system 
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described by p ps for any p ^ is a typical example. Then, 
can correlation between those local systems be regarded as 
classical one? An answer is found in the paper by Oppen- 
heim et al. (Ml (see also U) in which they made use of a 
class of states having a biproduct eigenbasis for a certain clas- 
sical/nonclassical separation. 

Definition 1. Let us consider a joint system consisting of sub- 
systems A and B with Hilbert space dimensions dffl and d\ B \ 
respectively. A complete orthonormal basis (CONB) consist- 
ing of eigenvectors, {\e.\ A ' B ^)}f = i dl \ of a density matrix of 
the system is a biproduct eigenbasis iff it is given by the direct 
product as {lef^)}, = {\ef)}f^ x {lejj? 1 }}^ where 

{\ e ^)}j an d {| e fc S ')}fe are eigenbases of individual subsys- 
tems. 

If a density matrix pi- 4 - 5 ! has a biproduct (BP) eigenbasis 

{|e^' B] >}i = {\ef)}gl x {\e[ B] )}i l :i it can be written in 
the form 

jk 

with coefficients < Cjk < 1 (Y2jk c i k = ^ state re P" 
resented by this density matrix is called a properly classically 
correlated state or, shortly, a classically correlated state J9[]. It 
is also called classical-classical state iflal . The state that can- 
not be represented by a density matrix in the above form is 
called a nonclassically correlated state. 

It is trivial to extend the above definition to a multipartite 
system. A density matrix having a (fully) product (FP) eigen- 
basis is written in the form: 

/',';; •"'* * if 1 c j ,..., x \ef)...\e^){ef\...(e^\, 

j,...,x=l,...,l 

(5) 

with the local CONBs {\ef ] )}fX {\e [ T ] )}tZ\ W 1 ' 1 is the 
dimension of the Hilbert space of a local system) and the coef- 
ficients < Cj a. < 1 Q3j x c j,—,x ~ 1) m use - This defi- 
nition separates the class of density matrices having a product 
eigenbasis (this is a non-convex set) from the class of density 
matrices having no product eigenbasis, as illustrated in Fig. 
[TJ The latter class is characterized by non-vanishing super- 
position (namely, non-vanishing off-diagonal elements) under 
local unitary transformations. 

In relation to the discussion on biproduct eigenbasis, Grois- 
man et al. Ill 111 recently introduced a measure of quantumness 
given by 

Pbp 

where the minimum is taken over all density matrices p^ B ' 
having a biproduct eigenbasis; F is any properly-defined 
distance function, such as the relative entropy function. 
Q(p^ A ' B ^) must be invariant under local unitary operations 
and must be zero for having a biproduct eigenbasis. 

They also suggested to use a special density matrix Pg^' 




FIG. 1 : Class of density matrices having a product eigenbasis (white 
regions in the class of separable density matrices) and the class of 
density matrices having no product eigenbasis (shaded region). 

(they called it Schmidt state) to define an easily-computable 
measure of quantumness F(p^ A ' B \ p^h^)- The density ma- 
trix Pg^h B ' i s created by keeping only diagonal elements of 
pl A , B \ under the special basis diagonalizing Tr^p^ 4,5 ' ® 

Tn/< , ,! - 

Indeed, it is a natural statement that a measure Al of non- 
classical correlation should satisfy the conditions: 

(i) M = if a system is described by a density matrix having 
a product eigenbasis (i.e., M = is a necessary but not suffi- 
cient condition for a state to have a product eigenbasis). 

(ii) M is invariant under local unitary operations. 

These conditions are considered to be prerequisite hereafter. 
In addition, one may test if a measure possesses either of (iii) 
full additivity, (iv) weak additivity, (v) subadditivity, etc. in 
the family of additivity properties. These properties are de- 
fined in the following way. They are desirable properties for 
measures of multipartite correlation and not exactly based on 
the additivity concept for bipartite correlation often seen for 
entanglement measures lfl2tl . Let us denote a measure of tri- 
partite nonclassical correlation by M m (er) where a is the den- 
sity matrix of an m-partite quantum system. First, the measure 
is fully additive if and only if 

M miXm2 (cri <g) cr 2 ) = M mi ((7i) + M m2 (a 2 ) 

with (Ti the density matrix of an mi-partite system and 02 the 
density matrix of an TO2 -partite system. Second, the measure 
possesses weak additivity if and only if 

M m »(^) = nl m (<7). 

Third, the measure possesses subadditivity if and only if 

M miXm2 (<Ti ®cr 2 ) < M mi (ai) +M m2 (a 2 ). 

In this paper, we introduce two measures of nonclassical 
correlation for a general multipartite system and numerically 
evaluate them for several examples. One of them, defined in 
SecHIl is similar to but different from the measure proposed 
by Groisman et al. and the other one, defined in Sec.[TTTJ is to- 
tally independent. This paper is organized as follows. In Sec. 
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Clare 



FIG. 2: System for which nonclassical correlation between Alice's 
part and Bob's part of a quantum state is discussed. A measure 
of nonclassical correlation is defined as a minimum uncertainty for 
Clare about the state after receiving their reports. 



HT1 we quantify a nonclassical correlation by a measure defined 
as the minimum uncertainty with respect to a joint system 
after we collect outcomes of particular local measurements. 
This measure satisfies the full additivity condition. The other 
measure will be introduced in Sec. [TTTJ which is defined in the 
following way: we consider the minimum distance between 
a genuine set and a mimic set of eigenvalues of a reduced 
density matrix of a local system on the basis of an artificial 
game in which one creates mimic eigenvalues of a reduced 
density matrix of a local system from eigenvalues of a den- 
sity matrix of a global system. The measure is defined by 
taking the maximum of this minimum distance over all local 
systems. It satisfies the subadditivity condition and a slightly 
stronger condition. We perform numerical computation of the 
two introduced measures for several examples in Sec.lIVIand 
compare them with negativity that is a common entanglement 
measure based on the separability paradigm. A discussion on 
definitions of nonclassical correlation is given in Sec.|V] Sec- 
tion [VI] summarizes the results of this paper. 

II. MEASURE OF NONCLASSICAL CORRELATION I 

We introduce the first of two measures of nonclassical cor- 
relation in this section. It is based on the paradigm in which a 
system described by a density matrix having a product eigen- 
basis is considered to possess only a classical correlation; in 
contrast, a system described by a density matrix having no 
product eigenbasis is considered to possess a nonclassical cor- 
relation. We have seen the form of a density matrix having a 
product eigenbasis for a bipartite system in Eq. and that 
for a multipartite system in Eq. ©. We will begin with the 
bipartite case. 



A. Bipartite case 

To quantify nonclassical correlation between distant sub- 
systems of a bipartite system, we consider the situation illus- 
trated in Fig. I2 Let us introduce Alice, Bob, and Clare. Alice 
and Bob have subsystems of a system and they are distant 
from each other. They can send reports to Clare. Alice/Bob 
can choose a complete orthonormal basis of her/his subsys- 
tem (basis {|e^ A ')}j for Alice and basis {lej. 3 ')}*; for Bob) 
for local (projective) measurements. Suppose that Alice and 



Bob use the observables Ma = Ylj j\ej )(e^ | and Mb — 

12k k\ e k ) ( e L I' respectively, and report the outcomes j and 
k to Clare. The probability that Clare receives j from Alice 
and k from Bob is Pjk = (e [ f \ {ef 1 \p^ A ^ \e [ f) \e[ B] ). The 
same process is performed for many copies of the same sys- 
tem shared by Alice and Bob without changing the bases ini- 
tially chosen. Then, minimum uncertainty (over all possible 
initial choices of local bases) that Clare has about pi- 4 '- 8 ! after 
receiving their reports is 

D(p^)= min ( -$> fe log 2 p jfe ) -Svn(p [A ' B1 ), 

local bases \ ' - / 

V > k / 

(6) 

where SVnCp'" 4 '" 8 ') is the von Neumann entropy. We em- 
ploy this quantity as a measure of nonclassical correlation. 
It is obvious that D(p^ A ' B ^) = for a density matrix with 
a biproduct eigenbasis. Otherwise, there is a possibility that 
D(p[^' B l) > 0. Thus a bipartite separable density matrix hav- 
ing no product eigenbasis (as well as a bipartite inseparable 
density matrix) possibly has a nonlocal correlation that may 
be quantified by D(p^ A ' B ^). A typical example is a density 
matrix shown in Eq. (O for < p < 1. For this density 
matrix, D(p^ A ' B ^) > holds. In addition, we should stress 
that D(p^ A ' B ^) is invariant under local unitary operations; this 
is clear from its definition in which we search over all local 
bases to obtain the minimum. 

One may find that the above process involving Alice, Bob, 
and Clare can be reconstructed in terms of CLOCC opera- 
tions HI |H] to define the same quantity: Consider the mini- 
mum discrepancy between the von Neumann entropy of the 
original state p^ 4 '- 8 ! and that of the state Clare can achieve 
by merging states received from Alice and Bob only after Al- 
ice and Bob locally use dephasing operations Ad{p^) = 

Si=i ( e i\p^ \ e i) \ e i) ( e i I where X is subsystem A or B, dJ- x ^ 
is the dimension of its Hilbert space, and {|e)i}i is the CONB 
of her/his choice. The minimum is taken over all choices 
of local CONBs. Then the discrepancy is equal to the mea- 
sure L>(p[^ s J). In this way, £>(p[ A,B l) can be related to 
the CLOCC protocol that was the base protocol for quantum 
deficit [9]. We should note that, in general, D(p^ A ' B ^) is not 
equal to quantum deficit; this is clear by comparing Eq. (11) 
of Ref. [9] with the above definition. The measure D(p^ A ' B ^) 
is equal to the quantum deficit in the case where the zero-way 
CLOCC protocol H] is considered |20ll . 

In addition, the measure D{p iA,B ^) is similar to but not in- 
cluded in the measure of quantumness defined by Groisman et 

al. HI 

in the sense that we search over all local bases while 
their measure is defined by taking a minimum over all classi- 
cal states (hence, over all possible combinations of eigenval- 
ues and local bases). 

B. Multipartite case 

It is straightforward to extend the definition Eq. (|6]l of D to 
that for general multipartite systems. Let us consider a den- 
sity matrix pi 1 '— > m l of an m-partite system. Consider local 
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CONBs {|e^)}j, . . . ,{|ei m ')} x . Then, a measure of nonclas- 
sical correlation is given by 



= , :n ! :i ~ V! /'/ - | "^/'/ - 

local bases \ ' * I 



(7) 



l,...,m]^ 



with 



Pi, 



[2] 



[1. 



An interpretation of this measure from an operational view- 
point is possible. Let L\ (L^) be the average code length of 
an optimal classical data compression, focusing on only di- 
agonal elements, acting on a register of qudits subsequent to 
local (global) unitary operations. Then the minimum discrep- 
ancy between L\ and L2 leads to the definition of the measure. 
The value of D(jfi''"' m ^) is zero if p^>~ > m l has a (fully) prod- 
uct eigenbasis. In addition, D(p^ L, '"' m ') is invariant under lo- 
cal unitary operations as is clear from its definition in which 
the minimum is taken over all local bases. Furthermore, we 
can prove that it is fully additive. 

The proof of full additivity uses only a general property of 
an entropy function in a particular form. Suppose that there is 
a system consisting of subsystems X and Y described by the 
reduced density matrices ax and ay, respectively. Consider 
an entropy function E of the product density matrix a x <3 
a Y and CONBs {|x)}, {\y)} of the reduced density matrices, 
written as 

E({\x)},{\y)},a x ®a Y ) 

= - (x\(y\a x <8> o Y \x)\y) \og 2 (x\(y\a x ®a Y \x)\y). 
\*)\v) 

We rewrite it as 

E({\x)},{\y)}, a x ® a Y ) = E({\x)}, a x ) + E({\y)}, a Y ) 



with entropies of subsystems 



J2\ x )(xWx\x) \og 2 (x\a x \ 



E({\x)},a x )-- 
E({\y)},<rY) = -E|„><»l°v|j/> ^og 2 (y\a Y \y). 

It is obvious that E({\x)}, {\y)}, ax ® cry) is minimized if 
and only if E({\x)}, ax) and E({\y)}, ay) are individually 
minimized. We conclude that D{a®r) = D(a) + D(r) holds 
for density matrices a and r accordingly. Thus the measure 
D satisfies the full additivity condition. 



C. Numerical method to estimate D 

A pure random search of local bases is a practical method 
to compute an estimated value of D(p'- 1 ''"' m ^), defined by Eq. 
(O, for a small size of a multipartite density matrix pi 1 '— > m l. 
Let us consider the bipartite case for clarity, where the density 



matrix is p 



[A,B] 



We employ an ad hoc random search algo- 



Algorithm: 

Continue to run the following trial: 
Trial: 

I. Generate a local CONB {|eJ A1 )} 3 X {lejjf 1 )}*. randomly by using the 
subroutine (written below) for both subsystems, A and B. 

II. Calculate the value of — J2jk Pj k 1°&2 Pj k ' 

A local CONB that minimizes — ^\ fc Pj k log 2 p^ k is chosen after many trials. 
The number of trials is typically 10 5 or 10 6 owing to a limitation of compu- 
tational resource. The measure D(p^ A ' B ^) is evaluated with the chosen local 
CONB. 

Subroutine (random generation of a CONB) — Gram-Schmidt Process : 

i. Generate a normalized (complex) d-dimensional random vector {v^ 
(here, d is the dimension of a reduced density matrix of the subsystem 
A or B for which we generate a basis). 

ii. For k = 2 to k = d, do 

1. Generate a (complex) ^-dimensional random vector \v k ). 

2. For I = 1 to I = k - 1, do 

\Vk) < \v k ) - ("l|«fc)|«l>- 

done 

3. Normalize 

4. if \Vk) is zero then goto i. 
done 

iii. Output { | v k )} d k=r 



FIG. 3: Algorithm to compute an estimated value of D(p' A,s ' 
the text for the symbols used herein. 



See 



rithm introduced in Fig.|3]in which a Gram-Schmidt process 



(see, e.g., page 108 of Ref. |QJ]) 

is utilized. In addition, it is 
needless to say that a similar algorithm can be used to estimate 
a value of D{p^ x '—' m ^) for a general m numerically. 

The random CONB generation is equivalent to a genera- 
tion of random unitary matrices acting on fixed CONBs {\j}} 
and {|fc)}. This is clear from the fact that for any local 
unitary operations Ua and Ub, we have Ua\j) — |e 3 -) and 
Us\k) — |ejfe). The probabilities pjfc are the diagonal elements 

of U\®U^ A ^U A ®U B . (It is now trivial to generalize this 
to the multipartite case.) Thus random matrix theories lfl4ll 
will be hopefully used to refine the algorithm. The present al- 
gorithm is still fast enough to estimate a value of D{p >1 '--- ,m ^) 
for a small number of qubits in a reasonable time. We show 
numerical results in Sec.lIVI 



III. MEASURE OF NONCLASSICAL CORRELATION II 

The previous measure, D, is defined as a minimum uncer- 
tainty about a global system after collecting outcomes of local 
measurements. Although the definition itself looks quite rea- 
sonable, we have to stress the difficulty to find the value of 
D because we need to try all possible product eigenbases to 
find the minimum, in principle. Strictly speaking, we need 
to try an infinite number of product eigenbases, but we rather 
find an estimate value by using a random search as we have 
seen. In this section, we will introduce the second measure 
derived from an artificial game to mimic a set of eigenvalues 
of a reduced density matrix of a local system by using a set 
of eigenvalues of a density matrix of a global system. The 
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FIG. 4: Illustration of the artificial game in which Alice (Bob) tries 
to find the eigenvalues of the reduced density matrix of her (his) part 
by the strategy described in the text. Clare knows the eigenvalues of 
the total system. 



advantage of this measure is that it can be calculated deter- 
ministically within a finite time. The disadvantage is its arti- 
ficial definition, which may be insignificant taking account of 
its advantage. 



may be given by 



F A (p^) = min 

partitionings 



2 C J 



(8) 

We may also consider the minimum uncertainty Fb (p^ A ' B ^ ) 
for Bob with respect to the set of eigenvalues for his local part 
in the same game. The larger one of their minimum uncertain- 
ties is then given by 



G(p [A > B] ) = m^{F A (p^),F B (p^)}. 



(9) 



This function can be used as a measure of nonclassical corre- 
lation since G(p^ A,B ^) = if p^ A ' B ^ has a product eigenbasis. 
In addition, G(p^ A ' B ^) is invariant under local unitary opera- 
tions according to the definition since local unitary operations 
preserve the eigenvalues of reduced density matrices of indi- 
vidual components and those of the density matrix of the total 
system. 



A. Bipartite case 

Consider an artificial game pertaining to Alice, Bob, and 
Clare. A system consisting of two parts, A (Alice's part) and 
B (Bob's part), is described by the density matrix p^ A ' B \ Al- 
ice and Bob cannot measure the system at all. Clare knows 
the set of all eigenvalues {ejk}fklfi ] where d\ A ^ and d\ B ^ are 
the dimension of the Hilbert space of Alice's part and that of 
Bob's part, respectively. Alice (Bob) wants to know the eigen- 
values of her (his) part. This setup is illustrated in Fig. [4] Let 
us concentrate on Alice's strategy. She asks Clare to send all 
the eigenvalues. Alice will partition x d\ B ^ eigenvalues 
received from Clare into d<- A ' sets 

To mimic an eigenvalue of the reduced density matrix of A, 
she sums up all elements of each set to make a set of 
mimic eigenvalues: 

{ei, • • • , e.j, ■ ■ ■ , e d \ A ]}, 

where each element is calculated as 

k=l 

The number of possible combinations to parti- 
tion d^ x d^ eigenvalues into d^ sets is 

d [A] x d [B\ \ / ^[A] d [B] \ / d [B] 

) { d [B] )■■■{ d [B] 

write the genuine eigenvalues of the reduced density matrix 
of A by {ej}j_l. The minimum uncertainty for Alice with 
respect to the set of these eigenvalues in the artificial game 



Let us 



B. Multipartite case 

An extension of the measure G to a multipartite case is 
straightforward. Let us consider an artificial game to find out 
eigenvalues of the reduced density matrix of a subpart from 
eigenvalues of the density matrix of the total system. Suppose 
that Kate has the fcth part of an m-partite quantum system. Let 
the dimension of the Hilbert space of the fcth part be tfM and 
that of the Hilbert space of the total system be dtot- She wants 
to know the eigenvalues {e^ }f=\ of the reduced density ma- 
trix of the fcth part. Kate receives d to t eigenvalues from Tony 
who knows the eigenvalues of the total system. Kate parti- 
tions them into d^ sets. Summing up elements in individual 
sets, she has d^ mimic eigenvalues {e^}. Thus a measure 
of nonclassical correlation for Kate can be 



F k (p 



') = min 

partitioning 



(ef log : 



* e 3 



log 2 ef) 



We may take the maximum over fc to have the measure 

G^pP'-'H) = maxF i; (p' 1 ' , "' m '). 

k 



(10) 



This is equal to zero if pi 1 ----'™] has a (fully) product eigenba- 
sis. In addition, it is invariant under local unitary operations as 
is clear from the fact that these operations preserve the eigen- 
values of a density matrix of the total system and those of 
the reduced density matrices of individual components. The 
measure does not satisfy the (full or weak) additivity condi- 
tion because the number of possible choices of partitioning 
eigenvalues grows rapidly as the system size grows. We can, 
however, prove its subadditivity. 

The proof for its subadditivity is accomplished in the fol- 
lowing way. Suppose the following story: there is a density 
matrix a of an m a -partite system having the set of eigen- 
values, {oia^Lj (here, d a is the dimension of the Hilbert 
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space of the system). With this set of eigenvalues, Kate has 
created a set of mimic eigenvalues of the reduced den- 
sity matrix of the fcth part. There is another density matrix 
t of an m T -partite system having the set of the eigenvalues, 
{AJb=r With this set of eigenvalues, Leo has created a set 
of mimic eigenvalues of the reduced density matrix of the 
kh part. Then a new game starts with a joint state a <B> r. 
The game for Kate/Leo is to make a set of mimic eigenvalues 
for her/his part by partitioning the product set of eigenvalues 
{«„} x {Pb} of a <g) r. Then, Kate (Leo) may make the same 
set of mimic eigenvalues as before because she (he) can make 
the set {a a } ({/?&}) firstly by partitioning {a a } x { /?(,}. Thus 
F k ((J <X> t) < Ffc(o-) and Fi(a ® r) < Fj(r) hold. Hence the 
next inequalities are satisfied. 

G{a <E>t)< max{G(cr), G(t)} < G{a) + G(r). 

In this way, the subadditivity G(a ® r) < G(er) + G(t) 
has been proved. We found that a slightly stronger condition 
G(cr <g> t) < max{G(<r), G(t)} is satisfied as a result. This 
property can be named submaximizability . 



IV. NUMERICAL RESULTS 



A. Examples for bipartite cases 



We compare the measures D and G with the negativity ifTM 



(here the map At is the 



Mmp [AB] ) = 

transposition map acting on B). 

The first example is the two-qubit pseudo-pure state 

p ps =p|VXVI + (i-p)i/4 

with \ip) = (| 00) + |ll))/>/2. We used a pure ( nume rical) 
random search of local bases, as introduced in Sec. Ill CI to es- 
timate a value of D(p ps ). The number of trials of local bases 
is 4.0 x 10 4 for each data point (namely, for each p in this 
example) of D(p ps ). (The number of trials is 4.0 x 10 4 for 
the other examples of two-qubit cases and is 4.0 x 10 5 for the 
examples involving an 8 x 8 density matrix.) Computation of 
a value of G(p ps ) is, in contrast, performed analytically. The 
eigenvalues of p ps are (1 + 3p) /4 and (1 — p)/4 with the mul- 
tiplicity three for the latter one. The eigenvalue of the reduced 
density matrix of a subpart is 1/2 with the multiplicity two. 
Thus we have 



1 -H 



l+p 



where 



H(x) = -a;log 2 2;- (1 - x) log 2 (l - a;) 



is the binary entropy function (0 < x < 1). 

Figure [5] shows the plots of D(p ps ), G(p ps ), and Af(p ps ) 
against p (0 < p < 1). The measures D and G reflect non- 
classical correlation for Vp except for p = 0. 
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FIG. 5: Plots of D(p pB ), G(p ps ), and Af(p pB ) against p. 



The next example is a mixture of Bell basis states, repre- 
sented by the density matrix 

p h =p\h)(bi\ + (l-p)\b 2 )(b 2 \ 

with = (|00) + |11»/V2 and \b 2 ) = (|01) + |10))/V2. 
This density matrix has a product eigenbasis if p = 0.5: 
p h (p = 0.5) = |(|+)|+)(+|(+| + |-)|-)<-K-|) with 

|±> = (|0>±|l»A/2. 

With the same basis search method as the previous exam- 
ple, we can estimate the values of D(p^). For this example, 
however, it can be analytically found: The von Neumann en- 
tropy of /?b is given by H (p) and the remaining term of Eq. (O 
becomes 1 irrespectively of the choice of a local basis. Thus 

D{p h ) = 1 - H{p). 

In addition, analytical computation of G{px>) is easy: the 
eigenvalues of p^ are p, l—p, and with the multiplicity two. 
The eigenvalue of the reduced density matrix of a subpart is 
1/2 with the multiplicity two. Thus, taking the minimum over 
partitionings of the eigenvalues of pb, we have 

G(p h ) = 1 - H(p). 

We find that D(p^) = G(pb) for this particular example. 

Figure [6] shows the plots of D(p^), G(ph), and Af(pb) 
againstp(0 < p < 1). All of these measures vanish at p = 0.5 
and have positive values for p ^ 0.5. Thus the difference 
among measures is not significant for pb- 

The third example is the density matrix (this is also of a 
2x2 system): 



1/2 -p 







p p 
p p 
1/2-p, 



(11) 



with < p < 1/2. This may be seen as a mixture of pure 
states 1 00), |11), and (|01) + |10))/\/2 with certain weights. It 
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FIG. 6: Plots of D(ph), G(pb), and Af(pb) against p. As we see in 
the text, D(pb) = G(pb) in this example. 
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FIG. 8: Plots of D(a[ A ' B] ), G{a [A ' B] ), and M{a[ A ' B] ) against b. 
Note that f^(a[ A ' B] ) is zero. 
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FIG. 7: Plots of D(a), G(a), and Af(a) against p. 



is separable for p < 1/4 because the eigenvalues of (7(g) At) a 
are 1/2 — 2p, p (with the multiplicity two), and 1/2. D(a) is, 
in contrast, non-zero unless p = as shown in Fig. [7] In the 
figure, we also plot 



min < 1 — H 



p , 1 - (2p) 



This is derived from the following values: the eigenvalues of 
a are 0, 1/2 — p (with the multiplicity two), and 2p; the eigen- 
value of Trscr and that of Tr^a are both 1/2 with the mul- 
tiplicity two. Interestingly, the shape of the curve of G(a) is 
similar to that of D(a). One drawback in using G{a) is that 
it vanishes at p = 0.25 although a has no product eigenbasis 
for this value. 

As the final example for the bipartite case, we consider 
the density matrix of a 2 x 4 system originally introduced by 



Horodecki |17]: 



1A,B] 



1 



76+1 



( b 














6 








\ 





b 














b 











b 














6 













b 




























1+6 
2 








Vl-b 2 
2 




b 














b 













b 














6 







Vo 





b 





Vl-6 2 
2 








1+6 
2 


/ 



(12) 



(0 < b < 1). This is known to be bound entangled for < 
b < 1, i.e., it is positive after partial transposition although it is 
inseparable for < b < 1. Hence Af(a\^' B ^) = 0. We find, in 
contrast, that D(al^' B ^ ) and G(a\^' B ^ ) are non-zero for b > 

as shown in Fig. [8] The measure G(cr[ A ' B ') is also plotted in 
the figure. This is accomplished by using the definition Eq. (0 
together with the following values. We have the eigenvalues 



of a 



[A,B] 



2b 



2b+l±V2b 2 -2b+l b 

146 + 2 ' 76+1' 76 + 1 

and (multiplicity three); 



(multiplicity two), 



the eigenvalues of Tr B (T^' B] : (36 + l)/(76 + 1) and 
46/(76 + 1); and the eigenvalues of Tv A a [ ^ B] : (36 + 1 ± 
Vl - 6 2 )/(146 + 2) and 26/(76 + 1) with the multiplicity 
two. 



B. Tripartite examples 

Numerical computation of D and G defined by Eqs. (0 
and ([Tol l, respectively, is easy also for a tripartite density 
matrix pi 1 * 2 ' 3 ! of three qubits. We will compare D(p^' 2 ' 3 ^), 
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FIG. 9: Plots of D(ppghz) and G(/9pghz) in comparison to those 
of the minimum and the maximum of the negativity for ppghz over 
all bipartite splittings (as functions of p). 
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FIG. 10: Plots of -D(cr[ 1,2 ' 3 '} 



and G(ct[ 1,2,3 ') in comparison to those 



of the minimum and the maximum of the negativity for a b 
all bipartite splittings (as functions of b). 



1,2,3] 



G^I 1 ' 2,3 !), and the minimum and maximum negativities of 
^[1,2,3] over a jj bipartite splittings. 

Consider the pseudo-GHZ (PGHZ) state: 

Ppghz = £>|^ghz)(^ghz| + (1 -p)l/8 

with IV'ghz) = (|000) + |111))/V2, the Greenberger-Horne- 
Zeilinger (GHZ) state. For this state, it is possible to find 
G(ppghz) analytically: The eigenvalues of ppghz are (1 + 
7p)/8 and (1 — p)/8 with the multiplicity seven for the lat- 
ter one; the reduced density matrix of any subpart has the 
eigenvalue 1/2 with the multiplicity two. Thus G(ppghz) = 
1 — H[(l + p)/2]. We can easily find that this equation holds 
for any number of qubits (larger than two) for the PGHZ state. 

Figure|9]shows the plots of -D(/0pghz), G(ppghz), and the 
minimum and maximum of the negativity for ppghz over all 
bipartite splittings against p (0 < p < 1). D and G reflect 
nonclassical correlation for Vp except for p = 0. 

The density matrix of the final example for a 2 x 4 bipartite 
case can be also interpreted as a density matrix of a 2 x 2 x 2 

[1 2 31 

tripartite case. Let us consider the density matrix a b ' of 
three qubits given by the matrix of Eq. d!21 l. Figure [lOl shows 
the plot of Z^crj, 1 ' 2 ' 31 ), the plot of G(cr[ 1,2 ' 3] ), and the plots 
of the maximum and minimum negativities over all bipartite 
splittings. The computation of G(cr| | 1 ' 2 ' 3 ') is performed by 

[1 2 31 

using the eigenvalues of <r l b ' ' ' (found in the previous sub- 
section) and the sets of eigenvalues of the reduced density 
matrices of subsystems k (k = 1,2, 3 and the following sets 
are given in this order): {(46)/ (76 + 1), (36 + l)/(76 + 1)}, 
{1/2, 1/2}, and {1/2, 1/2}. It is found that D(o-[ 1,2 ' 31 ) be- 
haves similarly to the maximum negativity over all bipartite 
splittings although the convergence values are different. The 
use of G(a| 1,2 ' 3 ') seems to be improper for this example be- 
cause its value is very small for 6 > 0.15 despite the fact that 

Fl 2 31 

a b : ' has no (fully) product eigenbasis for any 6. 



V. DISCUSSION 

We have considered a nonclassical correlation based on a 
commonly recognized paradigm other than the separability 
paradigm. We have introduced two measures of nonclassical 
correlation for the paradigm claiming that a multipartite sys- 
tem described by a density matrix having no product eigenba- 
sis possesses a nonclassical correlation. 

A motivation to consider a paradigm other than the separa- 
bility paradigm is connected to a conceptual investigation of 
the separability paradigm. Let us consider a bipartite system 
to simplify the discussion. The separability paradigm says 
that a system described by a bipartite separable density matrix 
can be prepared remotely when two distant persons (Alice and 
Bob) receive instructions from a common source. Hence there 
is no quantumness possessed by the system described by such 
a density matrix in the context of remote preparation. This 
assumes that a density matrix is a temporal average of instan- 
taneous density matrices or an ensemble average whose com- 
ponent density matrices are accessible independently because 
Alice and Bob prepare component states one by one. 

Nevertheless, this seems to be a non-mixing process be- 
cause Alice and Bob have accesses to individual instances or 
components; this is in contrast to usual processes in ensem- 
ble dynamics which are mixing. To make it a mixing process, 
they should lose their memories about time ordering of in- 
stances when the averaging is temporal averaging; in the case 
of ensemble averaging, they should lose their memories about 
indices of components. In contrast, Alice and Bob do not have 
any memory to lose when subparts of a system are distributed 
to them after a joint preparation of a quantum state. Thus, it is 
questionable to compare the amount and quality of correlation 
of a remotely-prepared state with those of a jointly-prepared 
state because the two contexts are different. 

The paradigm that we support in this work assumes that a 
density matrix having no product eigenbasis possesses non- 
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classical correlation. A discussion on the state preparation 
is not involved in its context. It is based on the problem as to 
whether or not off-diagonal elements of the density matrix of a 
multipartite system can be completely eliminated by local uni- 
tary operations. This paradigm has been widely recognized in 
the community while it is not considered to be a replacement 
of the separability paradigm. It is a highly conceptual prob- 
lem as to which protocol should be a base protocol to think 
about correlation. 

We have studied two measures D and G of non classical 
correlation based on the paradigm. The former is defined as a 
global uncertainty while the latter is defined as a local uncer- 
tainty. These two measures are quite different in their natu- 
ralness: the former can be interpreted as the minimum uncer- 
tainty for a global observer on the total system after receiving 
reports on local systems; the latter is based on some quite ar- 
tificial game and is, of course, unnatural. It is thus unexpected 
that plotted curves of these measures are sometimes similar 
to each other as we have seen in numerical results. We may 
enjoy the advantage to use the measure G, namely its easiness 
of computation. In contrast to D, which requires a numerical 
search for computing its value, G can be computed by consid- 
ering a finite number of eigenvalue partitionings. A drawback 
is that G is possibly small for a density matrix having many 
non-zero eigenvalues, as we have seen in Fig.[lO]for an exam- 
ple. It is easy to produce mimic eigenvalues of a subsystem 
close to genuine ones if the number of possible partitionings 
is large. This drawback is related to the fact that G is not 
(fully or weakly) additive but subadditive, in contrast to D 
that is fully additive. A measure with (full or weak) additiv- 
ity is more reliable to quantify nonclassical correlation as the 
system size grows. 

As we mentioned, there is a similarity between D and the 
measure defined by Groisman et al. lUlll . Let us consider a 
bipartite case. Suppose that we change the definition of D 
in the way that we choose the basis written as a product of 
two eigenbases of local systems A and B instead of search- 
ing the minimum over all product bases. Then, this redefined 
measure is included in the distance measures using the dis- 
crepancy between the given bipartite density matrix pt A ' B ^ 
and the Schmidt state p\£.^. The measure D is thus close 
to the measures defined as a discrepancy between p^ A ' B ^ and 
some specific density matrix having product eigenbasis (this 
can be the one that minimizes the distance, or some particular 
one to simplify the computing process). 

Finally, we take a brief look at an ongoing development 
in measures of nonclassical correlation. One of the mea- 
sures very recently proposed by Piani et al. [18] is designed 
to quantify nonclassical correlation in the same paradigm 



as presently employed. It is their measure Acc defined in 
a similar way as that of quantum discord [7]: Acc ls a 
discrepancy between the quantum mutual information lfl9ll . 
I(p^ A ' B ^), calculated for a bipartite density matrix p^ A ' B ^ and 
mmjvi A ^ B I[{M. a®Nb) p^' 8 ^] with two measurement maps 
M.A and Mb associated to POVMs. They also showed a 
straightforward extension of their theory to the multipartite 
case. An advantage of Acc 1S that it vanishes if and only 
if p^ A ' B ^ has a biproduct eigenbasis. A disadvantage is the 
difficulty of finding the minimum as is similar to the above- 
discussed disadvantage of D. 

There must be a large number of measures to quantify non- 
classical correlation for the paradigm; this is reminiscent of 
the dawn of entanglement measures. It is hoped that the 
paradigm will be studied extensively to extend another branch 
of quantum information science than the branch of the sepa- 
rability paradigm. 



VI. SUMMARY 

Two measures of nonclassical correlation have been intro- 
duced to support the paradigm claiming that a multipartite sys- 
tem described by a density matrix having no product eigen- 
basis possesses nonclassical correlation. The measure D has 
been defined as the minimum uncertainty about a joint system 
after we collect outcomes of particular local measurements. 
The measure G has been defined in the following way: con- 
sider the minimum distance between a set of mimic eigenval- 
ues and a set of genuine eigenvalues of a local system on the 
basis of an artificial game. The measure is defined by taking 
the maximum of this minimum distance over all local systems. 
We have shown that D is fully additive and G is subadditive. 
Numerical computations of D and G have been performed by 
using a random search of local bases and a non-probabilistic 
search of mimic eigenvalues, respectively. 
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